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PART  I 


DEFLECTION  OF  THICK  BEAMS  OF  MULTIMODULAR  MATERIALS 

C.W.  Bert"1"  and  F.  Gordaninejad* 

The  University  of  Oklahoma,  Norman,  Oklahoma,  USA  73019 


SUMMARY 


A  transfer-matrix  analysis  is  presented  for  determining  the  static  behavior 
of  thick  beams  of  "mul timodular  materials"  (i.e.,  materials  which  have  differ¬ 
ent  elastic  behavior  in  tension  and  compression,  with  nonlinear  stress-strain 
curves  approximated  as  piecewise  linear,  with  four  or  more  segments).  To 
validate  the  transfer-matrix  method  results,  a  closed-form  solution  is  also 
presented  for  cases  in  which  the  neutral -surface  location  is  constant  along 
the  beam  axis.  Numerical  results  for  axial  displacement,  transverse  deflec¬ 
tion,  bending  slope,  bending  moment,  transverse  shear,  axial  force,  and  loca¬ 
tion  of  neutral  surface  are  presented  for  multimodular  and  bimodular  models 
of  unidirectional  aramid  cord-rubber.  The  transfer-matri x  method  results 
agree  very  well  with  the  closed-form  solutions. 


INTRODUCTION 


In  1941  Timoshenko^  considered  the  flexural  stresses  in  bimodular  material, 
i.e.,  a  bilinear  material  having  different  moduli  in  tension  and  in  compres- 
sion.  Ambartsumyan”  in  1965,  introduced  the  terminology,  "bimodulus",  and 

For 

extended  the  concept  to  two-dimensional  materials.  Numerous  static  papers 

3 

appeared  after  this  work;  Marin  gave  the  effective  modulus  for  stiffness 
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of  bimodular  beam  undergoing  pure  bending. 

Small -deflection  bending  of  Bernoul li -Euler  beams  of  homogeneous,  bi¬ 
modular  material  was  treated  in  References  4-12.  Large  static  deflections  of 

1 1  C 

beams  of  bimodular  material  were  analyzed  in  References  13-14.  Kamiya 

considered  transverse-shear-deformation  effects  on  bimodular  beams  for  the 

1 6 

first  time.  Recently  Tran  and  Bert  treated  bending  of  thick  beams  of  bi¬ 
modular  materials  and  obtained  both  closed-form  and  transfer-matrix  solutions. 
To  the  best  of  the  present  authors’  knowledge,  no  previous  work  is  available 
in  the  context  of  multimodular  beams. 

MODELING  OF  THE  STRESS-STRAIN  CURVE 

Bert  and  Kumar^  recently  presented  experimental  stress-strain  curves  for 
unidirectional  cord-rubber  materials.  In  the  present  work  a  stress-strain 
curve  for  aramid-rubber  taken  from  [17]  has  been  linearly  approximated  by 
four  segments  (two  segments  in  tension  and  two  segments  in  compression).  For 
choosing  the  'break  points'  t  and  c  (see  Fig.  1)  the  area  between  two  fitting 
lines  and  the  experimental  curve  in  each  portion  has  been  minimized  (see 
Appendix  A).  To  find  comparable  moduli  for  the  bimodular  case,  one  has  to 
minimize  the  area  between  two  straight  lines  and  the  experimental  curve  (also 
see  Appendix  A) . 


THEORY  AND  FORMULATION 

Consider  a  rectangular-cross-section  beam  of  thickness  h  and  length  t  as 
shown  in  Fig.  2.  The  origin  of  the  Cartesian  coordinate  system  is  located 
on  the  mid-surface  of  the  beam  with  the  z-axis  being  measured  positive  down¬ 


ward. 
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1 .  Displacement  Field 

The  same  displacement  field  used  in  classical  Timoshenko  beam  theory  is 
implemented  here 

U(x,z)  =  u  ( x )  +  z;|,(x)  ,  W(x,z)  =  w(x)  (1) 


where  U  and  W  are  displacements  in  the  x  and  z  directions,  respectively, 
u  and  w  are  corresponding  displacements  at  the  midplane,  and  $  is  the 
bending  slope. 

2.  Stress  Field 


For  a  four-segment  approximation  of  the  normal  stress-strain  curve,  consider 

h  h 

ing  the  general  case  (i.e.,  when  -  ^  <  ac  ,  a£  <  j-),  the  following  stress 
field  has  been  considered  for  the  case  of  convex  bending  (see  Figs.  1  and  3) 
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r  r  t  t  C  t 

where  E i  ,  ,  E -j  ,  E2  »  G,  z^  ,  and  ^  are  material  constants,  ox  is  the 

axial  normal  stress,  s  is  the  axial  normal  strain,  y  is  the  transverse 

x  xz 

shear  strain,  :xz  is  the  transverse  shear  stress,  and  zn  is  the  location  of 
the  neutral  surface.  It  is  noted  that  this  material  is  linear  elastic  in 
shear.  Comparison  of  Figs.  1  and  3  leads  to 
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(6) 
(7) 


1  '  t  n 

efc  =  <(-h/2-  zn) 

=  <(h/2  -  zn) 


c  t 

where  £f  and  £f  are  the  final  values  attained  at  the  respective  compressive 
and  tensile  outer  fibers  and  <  is  the  curvature. 

Using  linear  strain  measure  and  the  strain  field  of  equations  (1),  one 
obtains 


E  =  U  =  U  +  Z’J  =  £  +  Z'l) 

X  ,  X  ,x  %x  ,x 


v  =  W  +  U  =  W  +  i!i 
’xz  ,x  ,z  ,x 


(8) 


Comparison  of  equations  (3)  and  (8)  gives 


u,x  '  -<2n  •  *,x  '  2  <9> 

Mote  that  (  )  denotes  d(  )/dx. 

j  X  * 

3.  Constitutive  Relation 

For  the  assumed  beam,  the  normal  and  transverse  shear  stress  resultants  and 
moment,  each  per  unit  width ,  are  defined  as 

h/2  h/2 

(N,0)  =  |  (ax,rxz)dz  ,  M  =  j  z^dz  (10) 

-h/2  ’  -h/2 


Using  the  assumed  stress  and  displacement  field  system  of  equations  (10) 
can  be  written  as  the  constitutive  relation  for  a  multi  modular  beam: 


i 
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where  A,  B,  D,  and  S  denote  the  respective  extensional,  flexural- 
extensional  coupling,  flexural,  and  transverse  shear  stiffnesses  defined  by 


,h/2 

( A , B , D )  =  j  (l,z,z2)E.(k)dz 
-h/2 

h/2 
2  f 

S  =  K  I  G  dz 

i 

-h/2 


1=1,2 

k=t,c 


(12) 


ABB  D 

Here,  the  stiffnesses  ,  and  are  not  present  in  linear  or 

bimodular  materials,  and  are  as  defined  in  Appendix  B.  In  equation  (12)  t 

and  c  denote  tensile-strain  and  compressive-strain  regions,  respecti vely . 

2 

The  quantity  K  is  a  shear  correction  coefficient  which  is  generally  taken 
to  be  5/6  for  static  loading  of  a  rectanqular-section  beam. 

4.  Equilibrium  Equations 

The  equilibrium  equations  for  transverse  distributed  loading  q(x)  can  be 
written  as 


N  =  0  ;  Q  +  q(x)  =  0 

,  A  >  A 


M,x  -  *  -  0 


(13) 


By  substitution  of  equation  (11)  into  equations  (13),  one  obtains  the 
following  equations  of  equilibrium  in  terms  of  the  generalized  displacements 


(A'u  -*■  B '  j  )  =  0 

,x  ,x\x 

[S(w  +  v)]  =  -q(x) 

♦  A  ,  A 


(B"u  +  D'-:  )  -  S(w  +  v)  =  0 

,A  1  A  J  A  )  A 


(14) 


where 
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♦  c* 

B1  ■  B  +  C8 
B"  =  B  +  CB 


(15) 


CLOSED-FORM  SOLUTION 


A  closed-form  solution  can  be  obtained  only  when  the  stiffnesses  and  thus 
neutral -surface  position  (z  )  do  not  depend  on  x.  Therefore,  neutral- 
surface  location  [18]  must  be  constant 

2n  =  'u  x  =  constant  (15) 

Using  equation  (11)  one  is  able  to  express  u  and  >;  in  terms  of  N,  M,  and 

>  X  5  X 

the  stiffnesses  as  follows: 
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Combining  equations  (16)  and  (17),  one  obtains 

zn  =  (B'M  -  D ' N ) / ( A ' M  -  B"N)  (18) 


It  is  obvious  that  z  =  const,  when  N  =  0.  Thus,  for  these  special  cases 


(see  Appendix  C) 


zp  =  B'/A'  =  constant 


(19) 


Now,  equilibrium  equations  (14)  can  be  simplified  as  follows  [19]: 
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,xx  -  S(w,x  +  ■>  =  0 


(20) 
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The  genera]  solution  for  equations  (20)  can  be  written  as  follows: 
u ( x )  =  d1  +  d2x  +  ~  C4x2  +  u  (x) 

*(x)  =  -C2  +  6(--,^-rA-—  C4 -  2C3x  -  3C4x2  +  u»p(x)  (21) 

w(x)  =  C1  4  C2x  +  +  C4X3  +  wpu) 

where  u  ,  ,  w  are  particular  solutions  (see  Appendix  D)  and  ,  C2>  C^, 

C4,  d^ ,  and  d2  are  arbitrary  constants  determined  by  the  boundary  conditions 
of  the  beam.  The  following  boundary  conditions  have  been  considered  for 
closed-form  solutions: 

1.  Hinged-Hinged  (free  to  move  axially  at  x  =  L) 

u(0)  =  N(i)  =  0  ;  M(0)  =  M(i)  =  0  ;  w(0)  =  w(t)  =  0 

2.  Clamped-Free 

u(0)  =  N(0  =  0  ;  *(0)  =  M(i)  =  0  ;  w(0)  =  Q(0  =  0 

3.  Clamped-Clamped  (free  to  move  axially  at  x=L) 

u  ( 0 )  =  N  (i)  =  0  ;  *(0)  =  *U)  =  0  ;  w(0)  =  w  (l)  =  0 

The  values  of  constants  C-|,  C2>  Cj,  C4>  d-j ,  and  d2  are  listed  in  Appendix  D. 

TRANSFER-MATRIX  SOLUTION 

As  it  has  been  shown  in  [20,21],  in  the  transfer-matrix  approach,  the  beam 
is  divided  into  Ns  elements,  each  of  which  is  assumed  to  be  of  mass  m  and  con¬ 
centrated  at  the  center  of  mass  of  the  element.  The  mass  center  of  each 
element  is  called  the  station.  The  stations  are  separated  by  fields  which 
are  taken  to  be  massless  and  contain  all  of  the  stiffnesses  of  the  beam. 

At  the  end  points  of  the  beam  there  are  two  half  fields  of  length  n',/2  (see 


8 


Fig.  4),  and  between  these  half  fields  there  are  N^  stations  separated  by 
(Ns  -  1 )  full  fields  of  length,  m,  where  At  =  ;/N  and  i  is  the  length  of 
the  beam.  By  writing  the  equilibrium  equations  for  each  station  and  each 
element  and  connecting  the  elements  by  transfer  matrices,  one  transfers  the 
generalized  displacements  (u,w,y)  and  the  forces  (N,Q,M)  from  the  left  side 
of  the  beam  to  the  right  side. 

Since  the  same  procedure  used  in  [16]  has  been  used  here  to  derive 
transfer  matrices  and  state  vectors  (u,w,q,N,Q,M)"'\  the  readers  are  referred 
to  this  reference.  Note  that  since  the  present  work  deals  with  multimodular 
material,  some  changes  in  the  field  matrix  are  necessary  (see  Appendix  E). 

In  the  calculation  of  the  stiffnesses  for  the  cases  where  the  axial  force  is 
not  zero,  the  neutral -surface  locations  and  the  corresponding  distances  to  the 
"break  points"  in  the  cx  vs  z  curve  (ac  and  at)  are  not  constant  and  not 
known  a-priori.  Therefore,  an  iterative  technique  has  been  employed  to  com¬ 
pute  the  neutral -surface  locations  also  ac  and  at>  One  must  first  assume 
(2N  +2)  sets  of  values  of  zn,  a  ,  and  at  and  then  compute  the  stiffnesses 
and  solve  the  governing  equations  for  the  state  vector.  Finally,  by  using 
equations  (18),  (C.3),  and  (C.4),  compute  new  values  of  z  ,  ac»and  a^. 
Obviously,  if  the  assumed  and  computed  sets  of  z  ,  ac>  and  at  are  in 
sufficiently  close  agreement,  the  problem  is  solved;  otherwise,  assume  the 
calculated  set  zn,  ac>  and  at  and  repeat  the  procedure. 

NUMERICAL  RESULTS 

In  the  following,  numerical  results  are  presented  for  a  thick  beam  with 
a  rectangular  cross  section  and  constructed  of  multimodular  material  (see 
Table  1):  unidirectional  aramid  cord-rubber,  which  is  used  in  the  tire 
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industry.  Various  boundary  conditions  and  loading  conditions  were  investi¬ 
gated  (see  Tables  2  and  3).  In  the  transfer-matrix  analysis,  twenty-five 
elements  were  used.  Each  element  was  of  length  0.32  in.  for  dimensional 
cases  and  dimensionless  length  of  0.04  for  nondimensional  cases.  The  shear 
correction  coefficient  was  taken  to  be  5/6. 

For  all  cases  considered,  the  computations  are  carried  out  for  axial 
elongation  u  (or  u  =  u  E^/q^),  transverse  deflection  W  ( or  W  =  W  E^/q^O, 

bending  slope-,  (or  ;  =  .  ^^o^  axia^  f°rce  N  (or  N  =  N/q  :),  shear  force 

Q  (or  Q  =  Q/qQ0,  bending  moment  M  (or  M  =  M/qoc:)  and  neutral -surface  loca¬ 
tion  zn  (or  zn  =  zn /h ) ,  where  u,  W,  N,  Q,  M,  and  zn  are  nondimensional 

parameters . 

Due  to  lack  of  comparable  results  in  the  literature,  comparisons  are 
made  between  the  closed-form  solution  (CFS)  and  the  transfer-matrix  solution 
(TMS)  developed  here.  Excellent  agreement  between  CFS  and  TMS  for  the 
twenty-five  element  model  has  been  achieved  and  still  it  can  be  improved  by 
increasing  the  number  of  elements.  For  most  of  the  results,  the  error  is 
less  than  2%.  Figure  5  contains  the  plots  of  the  dimensionless  transverse 
deflection  (W)  versus  dimensionless  position  (X  =  x/z)  for  Case  11  for  multi- 
modular,  bimcdular,  unimodular,  and  average-modular  [i.e.,  E  =  (E^  +  E^  )/2] 

cases,  where  ;/h  =  10. 

As  one  can  see,  there  is  a  considerable  difference  between  transverse 

defection  of  multimodular  and  oimodular  models  on  one  hand  and  unimodular 

and  average  modular  models  on  the  other.  Note  that  Case  11  is  a  special 

case  because  both  ends  are  not  free  to  move  and  one  expects  axial  force  to 

be  developed  due  to  bending-stretching  coupling  caused  by  bimodular  action. 

However,  the  commuted  axial  force  is  close  to  zero  which  means  z  is  constant. 

n 
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To  validate  TMS,  in  Figs.  6  through  9,  a  comparison  is  made  between 
TMS  and  CFS.  Behavior  of  Cases  1,  4,  10,  and  11  was  studied  for  the  multi- 
modular  model  considering  different  dimensionless  parameters  ( n/h  =  5,  10,  15, 

M]  =  -  1.0,  ^  =  -  1.0,  ^  =  -  1.0). 

In  Table  4,  for  a  specific  beam  (see  Table  1),  dimensioned  comparisons 
have  been  made  between  multimodular  and  bimodular  models.  Tables  5,  6,  7 
and  9  again  show  the  validity  of  TMS  while  they  present  the  computed  y'esu 
for  Cases  3,  5,  6,  7,  and  8. 

Since  closed-form  solutions  are  not  available  for  the  complicated 
boundary  conditions  considered  in  Cases  9,  10,  and  11,  only  transfer-matrix 
results  are  presented  for  these  cases.  In  Figs.  10-13  the  behavior  of  a 
clamped-free  beam  (with  applied  moment  and  axial  and  shear  forces  at  the  free 
end)  and  a  clamped-cl amped  beam  under  a  uniform  load  is  investigated.  See  also 
Table  10.  It  is  of  particular  interest  to  note  that  in  Fig.  10,  the  sign  of 
the  deflection  depends  upon  the  t/h  ratio,  the  crossover  point  being  at  t/h  =11. 

For  cases  where  axial  force  is  zero,  the  neutral -surface  location  is 
constant;  otherwise  it  varies  along  the  beam  length.  In  Figs.  14-17  the 
shapes  of  neutral -surface  curves  for  multimodular  and  bimodular  models  have 
beeh  shown  for  several  cases. 

CONCLUSIONS 

Analyses  of  the  bending  deflection  of  multimodular  thick  beams  with 
rectangular  cross  section  based  on  shear-deformable-beam  theory  are  presented. 

In  this  study,  both  dimensionless  and  dimensioned  results  of  transfer-matrix, 
as  well  as  closed-form,  solutions  for  a  rectangular  multimodular  beam  of 
aramid-cord  rubber  are  presented.  The  transfer-matrix  and  the  closed-form 
solutions  are  found  to  agree  very  well. 
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Results  of  analysis  of  bimodular  and  multinodular  models  show  that 
there  is  not  a  drastic  difference  between  the  two  models.  Although  the 
mul timodular  model  is  a  better  one  for  approximating  the  stress-strain 
curve,  the  bimodular  approximation  is  less  complicated.  Closed -form 
solutions  are  available  only  for  a  number  of  loading/boundary  conditions 
(in  which  the  axial  force  is  identically  zero),  but  the  transfer-matrix 
method  can  be  applied  to  more  complicated  geometry,  loading,  an!  boundary 
conditions.  In  this  work,  results  for  several  boundary  and  loading  condi¬ 
tions  are  investigated.  The  transfer-matrix  method  is  found  to  be  very 
effective  in  terms  of  computational  time  and  also  gives  results  which  agree 
quite  well  with  the  closed-form  solutions. 
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APPENDIX  A 

FITTING  MINIMIZED  CURVES  TO  THE  STRESS-STRAIN  CURVE 


1 .  Multimodular  Case 

Consider  the  nonlinear  stress-strain  curve  shown  in  Fig.  A.1.  For  any  arbi¬ 
trary  point  (£t,ot)  in  the  tension  region  (e  >_  C),  there  are  two  straight 
lines  such  that 


<j(e) 


f (at/£t)e 

a^t)/(£t  -  £^)](e  -  e^)  + 


(A.1) 


The  equation  of  a  stress-strain  curve  as  expressed  in  [17]  is 

c( c)  *  Kcn  ;  £  >  0  (A. 2) 

where  K  and  n  are  constants  depending  on  the  material.  To  find  the  proper 
"break  point"  the  area  between  the  approximated  curve  g(c)  and  the 

actual  experimental  curve  j(c)  has  to  be  minimized.  The  mentioned  area  can 
be  expressed 


A  = 


t 

c 

j  [a,  (c)  -  o(e)]de j  + 

i  1  I 

0 


t 

['f 

!  [g2(E)-:(r)]de 

J  L 

t 


(A. 3) 


Substitution  of  equations  (A.1)  and  (A. 2)  into  equation  (A. 3)  and  taking  the 
integrations  gives 


A  =  jl  -t, 


K 

n+1 


,  n+1  i 

(cl)  ! 


,t  +  ot)(: 


n+1 


-  (cl) 


n+1 


(A. A) 

By  searching  in  the  region  of  =  ( 0 , £ t ’ )  x  [0,-^),  one  is  able  to  find  a 
point  (£t,c^)  such  that  A  is  minimized  locally.  Note  that  a  few  other 
methods  (e.g.,  least-squares  method)  have  been  tried  but  it  turned  out  that 


the  absolute  minimum  point  was  outside  of  the  region  2. 

2.  Bimodular  Case 

For  this  case,  the  least-squares  method  has  been  used.  As  shown  in  Fig.  A. 2, 


there  is  a  line  such  that 


I 


Ef 


(A. 5) 


can  be  minimized  in  Here,  E^  is  the  slope  of  that  line.  By  taking  the 
derivative  of  equation  (A. 5)  and  equating  it  to  zero,  one  has 


By  solving  equation  (A. 6)  for  ,  one  obtains 

n-1 


t  =  3K 

3  n+7 


Uf1) 


(A. 6) 


(A. 7) 


For  example,  for  arami d-rubber  in  the  tension  region  (see  [17]),  the 
following  parameters  are  found: 


nt  -  1.22 

Kt  =  1 . 1  x  106  psi 

-;ft  =  0.029 

r  t  _  (3)  (1 . 1  x  TO6) 
bb  1.22'+  2 


(0.029) 


1.22-1 


0.47  x  106  psi 


An  analogous  calculation  can  be  applied  for  the  compression  side  of  the 
bend,  i . e . ,  £bC  can  be  found,  provided  that  Kc,  nc,  and  zfc  are  known. 


[ 
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APPENDIX  B 

THE  BEAM  STIFFNESSES  FOR  RECTANGULAR-SECTION 
BEAMS  OF  MULTIMODULUAR  MATERIALS 

For  the  assumed  four-segment  model,  there  are  two  different  bending  cases 
in  general,  convex  downward  and  concave  downward  bendina.  In  convex  down¬ 
ward  bending,  the  top  layer  of  a  beam  is  in  compression  and  the  bottom  layer 
in  tension.  Conversely,  in  concave  downward  bending,  the  top  layer  of  the 
beam  is  in  tension  and  the  bottom  layer  is  in  compression. 

Depending  on  the  location  of  zn,  a  ,  and  a^  in  cx  vs  z,  eight  different 
cases  might  occur.  For  example,  for  convex  downward  bending,  consider  the 
case  when  zn,  a  ,  and  at  are  in  the  range  of  -h/2  and  h/2  (see  Fig.  3). 
Substitution  of  equation  (2)  into  equation  (10)  and  using  equations  (3),  (4), 
and  (5)  leads  to 


•zn  rat 

[<E i c ( a  -  z  )  +  <E9c(z-a  )]dz  +  I  <E,c(z  -  zjdz  +  f  <E,t(z-z)dz 


-h/2 


1  v  n'  J  1 


+  j  [<E1t(at  -  zn)  +  <E2t(z  -  at)]dz 


'  c  r  r  r  n  r 

M  =  [<E1  (ac  -  zn)  +  <E2  (z  -  ac)]z  dz+  !  <E]  (z  -  zp)z  dz 

-h/2 

,at  h/2 

+  |  <E1t(z-zn)z  dz  j  [<E-,t(at  -  zn)  +  kE^U  -  ) ] z  dz 


( B .  2 ) 
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Equations  (B.l)  and  (B.2)  can  be  written  in  the  following  form 


A  .  fZn  ,at 

(-<zn)  {[,  e2  dz  +  i  e,  dz  +  j  eA 


2  dz  +  j  E]  dz  +  j  dz  +  j  E2  dz] 


,ac  rac  h/2  .h/2 

C-  |  E2zdz  +  j  E-|Cdz  +  ;  E-|ldz  -  f  E^dz]} 

‘h/2  -h/2  a.  a. 


,ac  ,zn  _at  h/2 

+  (<)  {[;  E?cz  dz  +  |  E-j  Cz  dz  +  j  E^z  dz+  ^  E^ 

-h/2  a  zn  a. 

c  n  t 


z  dz] 


rac  c  <-ac 

+  n  E2  acdz  -  I  E1 


h/2  h/2 

ElCacdz  +  j  E]tatdz  -  |  E2tatdz]} 


,ac  ,zn  at  h/2 

(~<zn)  {[j  E2Czdz+!  E  i c  z  dz  +  J  E^zdz+f  E2tz  dz] 

-h/2  ac  zn 

•ac  r  rac  A2  fh/2 

+  [-  j  E2cz  dz  +  |  E^zdz+j  E^zdz-j  E2tz  dz] 


•ac  -zn  -at  -h/2 

+  (-•)  {[j  E2cz2dz  +  ;  E1cz2dz  -  j  E]tz2dz  +  j  EgVdz] 

-h/2  a_  z 


,ac  ,ac  .h/2  _h/2 

+  [j  E2°acdz  -  j  ElCacdz  +  j  Eltatdz  *  I  E2tatdz]} 


Combining  equations  (9)  and  (11),  one  gets 
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+  <B  1 

(B.5) 

■  <-'2n>B" 

+  <D' 

(B.6) 

Comparison  of  equations  (B.3)  and  (B.5)  with  equations  (B.4)  and  (B.6)  and 
considering  equations  (12)  and  (15),  one  finds  that 


h/2 


CAN  =  j  (E^-E^dz. 
-h/2 


(E1t-E2t)dz 


:bN  -  J  (E^-E^dz  + 

-h/2 


h/2 


(E1t-E2t)atdz 


ac  h/2 

=  |  (E1C-E2c)z  dz+  j  (E1t-E2t)zdz 


-h/2 


M  _  f  c  c 


h/2 


(E,'  -  E ,  ) a _z  dz  + 


'D  J  v  “2  "1  'V 

-h/2 


(E] 1  -  E2t)atz  dz 


(B.7) 


As  mentioned  before,  eight  cases  may  occur  depending  on  the  location  of  zn, 
ac>  and  a^.  These  cases  have  been  analyzed  as  the  same  as  the  general  case  as 
follows  (for  convex  downward  bending) 

Case  1 : 


A1  *  hE1 1 

S'  ■  ^(E/.E^) 

B"  =  0 

0'  =  h3E2V24 


(B.8) 


(B.  9) 
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Case  2: 


E  l-  1  +  E  t(, 
4  '1  2  Ux 


E 

v  1  "x 


at  <  z  <  h/2 
-h/2  <_  z  <_  at 


(B. 10) 


A1  =  hE-j 1 

B*  =  -  (E]t-E2t)(h/2-at)2/2 
B"  =  0 

D 1  =  [h3(E1t  +  E2t)/S  +  at3(E2t- E] t)/2]/3  +  h2afc(E1 1  -  E2t)/8 


(B.  1 1 ) 


Case  3: 


^  p  t _  t  +  p  t /  t\ 

E1  “1  4  1  X  M  ] 

r  t 

E1  ex 

E  c  ■ 

.  1  x 


at  1  z  1 

z„  <  z  <  a. 
n  -  -  t 

-h/2  <  z  <  2^ 


(B.ll) 


A'  =  h(E1C+E1t)/2  +  zn(E1c  -  E1t) 

B1  =  [h2(E2t  -  E]c)/2  +  at(h-at)(E1t-E2t)  +  zp2  (E,c  -  E^JJ/2 
B"  =  (E^  -  E1c)(h2/4  -  rn2)/2  (B.  12) 

O'  =  [n3(E1c  +  E2t)/8  +  at3(E2t-E]t)/2  +  -  E]  t)]/3 

+  h2at(E1t  -  E2t)/8 


Case  A: 

This  case  is  the  general  case  (see  Fig.  3)  which  has  been  discussed  in 
detail  earli_.  in  this  Apppndix. 

Case  5: 


Vs 


-  J 


1  x 

c  _  C  +  p  c2  . 

1  '  1  b2  v'-x  "'1  } 


zn  -  Z  -  h/2 

ac  I  z  i  zn 
-h/2  <  z  <  ac 


(B.  13) 
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A1  =  hCE^+E^)  +  zf|(E1c-E1t) 

B'  =  [h2(E1t  -  E2c)/4  +  zn2(^c-t*)  +  ac(ac+h)(E1c-E2c)]/2 
B"  -  (E^  -  E1C)(h2/4-zn2)/2  (B.  14) 

0’  =  [h3(E1t  +  E2c)  +  2n3(E1C  -  E/)  +  ac3(E1c-E2c)/2]/3 
-  h3ac(E1c-  E2c)/8 

Case  6: 

r  c  C 

E1  £x  ac  -  2  -  h/2 

cx  H  ‘  (B. 15) 

,EiCmC  +  E2C(ex-^C)  -h/2  <  z  <  ac 


A'  =  hE-j C 

B'  =  (E-,c  -  E2C) (h/2  +  ac)2/2 
B"  =  0 

D'  =  [h3(E1C  +  E2c)/8  +  ac3(E1C-E2C)]/3  -  h^E^  -  E2C)/8 

Case  7: 

ax  3  E1C  1°  +  E2C^£x  '  e]C)  -h/2  1  z  1  h/2 

A'  =  hE-|C 

B'  -  hac(Elc-E2c) 

B"  =  0 

0'  =  h3E2C/24 


Case  8: 


l _ 


x 


E 


zn  i  z  i  h/2 
-h/2  <  z  1  zn 


(B.  17) 


(B.  18) 


(B. 19) 
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A'  =  h(E]t+E1C)  +  zn(E1C  -  E1t) 

B'  =  (E^  -  EiC)(h2/4  -  zn2)/2 

B"  =  (E^-  EjC)(h2/4  -  zn2)/ 2 

D1  =  [h3(E-|t  +  E-|C)/8  +  zn3(E1C  -  E1 1)]/ 3 


For  concave  downward  bending,  one  is  able  to  derive  similar  equations  for 
stiffnesses  by  converting  as  follows: 


E 

E 


2 

1 


c 

c 


a 

a 


t 

c 


APPENDIX  C 

COMPUTATION  OF  z  ,  a  ,  AND  at 

For  mu  1 timodnular  beams,  the  following  equation  is  not 
mine  the  neutral -surface  location  zR 

B'M  -  D'N 
2n  ~  A'M  -  B"N 


sufficient  to  deter- 


(C.l) 


even  for  cases  where  N  =  0 

zn  =  8  7A'  ( C .  2 ) 

Two  more  equations  are  needed  for  computing  zn  because  the  stiffnesses  are 
not  only  dependent  on  zn  but  they  are  functions  of  ac  and  as  well. 
Dividing  equation  (4)  by  equation  (6)  and  equation  (5)  by  equation  (7)  3nd 
solving  for  ac  and  a^,  one  can  get  (for  the  convex  downward  case) 


■  (-,c/.fc)  (h n  *  in)  -  zn 

(C .  3 ) 

-  (-:,*/■:/)  (h/2  -  Z„)  -  Z„ 

(C.4) 
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For  the  concave  downward  case 


ac  ■  (c,c/efc)  (h/2  -  z„>  *  zn 
at  ’  C/2  *  z„>  '  z„ 


The  system  of  nonlinear  equations  (C.l),  (C.3),  and  (C.4)  for  the  convex 
downward  case,  or  equations  (C.l),  (C.5),  and  (C.6)  for  the  concave  down¬ 
ward  case, can  be  sol ved  by  using  iteration  of  the  Gauss-Seidel  type  [22]. 

APPENDIX  D 

ARBITRARY  CONSTANTS  AND  PARTICULAR  SOLUTIONS 
The  values  of  constants  C-| ,  C2,  C^,  C^,  d-j,  and  d^  for  the  various  boundary 
conditions  considered  are  listed  below. 


Hinged-Hinged  (free  to  move' axially  at  x=L) 


C,  =  -  up(0) 

C2  =  -  (C3;  +  C4;.2)  -  [wp( ;)  -  wp(0)]/: 


c,  =  [vn  y(i)/2]  -  3C. L 


C4  =  6T(A’D‘  -  B 1  S " )  B'fUp,x^)  '  Up,x^0^  +  D'f-p,x(-) 


d1  -  -  up(0) 


W°)] 


d,  =  -  [un  VU)  +  ( SB ’ /A 1 ) C ,  t ] 


2.  Clamped-Free 

C1  =  -  wp(0) 

C2  =  -  wp,x(0  *  V')  +  -P(0) 

C3  =  vp,x^/2  +  2  ( B '  B" A 1 D 1  )  ^wp,x('^  + 


CA  =  6 ( B  ' 3"  -  A 1 D '  )  'C2  '  ;p(G'n 
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J1 


u  (0) 
P 


B'Se 


d2  "  ’  Up,x(;)  +  B  'B"  -  A'D'  ^wp,x(^  +  vp('^ 


3.  Clamped-Clamped  (  free  to  move  axially  at  x=  L) 
C,  =  -  w  (0) 

1  pv  ' 

6(B'B"  -  A ' D '  ) 


SA* 


C4  +  •p(°) 


(D.3) 


C3  =  [un ( 0)  -  u„(0  -  (3B7A‘)C nr}/2- 


C,  = 


P  P 

SA' 


4  "  SA' 7  -  lTTB'B"  -  A '  D ' )  "p 


rnTT  ^ A  ( 0 )  +  vn(0  -2[w  (0)  -w 


dl  =  -  up(0) 

d2  =  [ u D ( 0 )  -  u  (0  -  (3B ' /A '  )C4 ■'“]/>: 


The  particular  solutions  for  uniform  and  sinusoidal  normal  lead  are  as 
1 i sted  below. 

For  uniform  normal  load  q(x)  =  qQ: 

5 '  o 

t  \  -  o  3 

V*  '  STA'O'  -  3' 3^7 

A'q. 


.  (x)  =  -  x  - 
P  i 

A'q 

/  \  0 
V*  =  24(  A '  D '  -  B'FT  X 


STAG'D B  r‘g^TT 


:or  normal  load  q(x)  =  q  sin  -x,  where  *  =  n-/.': 
3 '  q 

=  cos  ,X 


= 


A'q, 


T^A'D’  -3'8'T 


-x  POS  :X 


n  qO  rl  A' 

Vx  '  s  xTA'd' -B’B'T  s,n  lXJ 


(0.4) 


(0.5) 
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For  normal  load  q(x)  =  qQ  cos  ^x: 


B'q0 

Up(x)  ( A ‘ D ‘  -  S  1 8" )  Sin  aX 

A 1  q 

■Pp(x)  =  aj(A.D.  sin 


,  ,  _  qO  rl  X  A1  1 

V*  '  '-S  ar  ( A 1 D 1  -  B  1 8" )  C0S  aX 


(D.6) 


APPENDIX  E 
TRANSFER  MATRICES 

The  equilibrium  equations  for  each  station  can  be  written  in  matrix 
notation  as  follows 


f u  ^ 

R 

n 

0 

0 

0 

0 

0 

0 

f 

w 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

! N 

►  . 

0 

0 

0 

1 

0 

G 

0 

< 

0 

0 

0 

0 

0 

1 

0 

% 

M 

0 

0 

0 

0 

0 

0 

1 

l'. 

0 

0 

0 

0 

0 

0 

1 

— J 

(E.l) 


Q 

m  : 


where  qs  is  the  concentrated  load  at  each  station, 
equation  (E.l)  is 

[S],  =  [Ts].[S]. 


In  more  compact  form 


(E.2) 


The  matrix  [T<-]  is  known  as  the  station  matrix.  In  matrix  notation  the 
i 

equilibrium  equation  for  each  field  under  a  distributed  load  q(x)  is 
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j 

0 


Values  of  '<  and  K  for  various  loadings  are  listed  in  Table  2.  Equation 
m  q  3 

(E.t)  also  can  be  written  as 

[S]L  =  [T-]  [S]R 
i  1  J  i  i 

Tne  matrix  [T^].  is  called  the  field  matrix. 


(E.5) 


Parameters  Elastic  ^opeKies 
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Table  1.  Elastic  Properties  and  Geometric  Parameters 
for  an  Aramid-Cord  Rubber  Beam 


Longitudinal  Young's  Modulus, 
psi  x  10-6 

Longi tudinal -Thi ckness  Shear 
Modulus,  psi  x  10"~ 

Tension 

Compression 

M 

mtm 

E  C  E  c 

tl  t2 

G 

0.580  0.420 

0.032  0.010 

0.537 

B 

C 

0.470 

0.180 

0.537 

r~ 

c. 

E 

0.275 

0.537 

0.244 

0.537 

Beam  length  8.0  in. 
Beam  depth  (thickness)  0.6  in. 
Beam  width  1.0  in. 


I _ 1 - - - J 

M  Mul timodular ,  B  n-  Bimodular,  U  Unimodular,  A  %  Average  Modular 


'able  2 

1 - 


29 


Summary  of  Casas  Considered 


Case  Ccj nda.-y  Condition  anc 

No.  Load  Position 


Sco.odary  Condition  an 
Load  Position 
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Table  3.  Values  of  K  and  K  for  Various  Loadings 

m  q  3 


Table  5.  Closed-Form  and  Transfer-Matrix  Solutions 
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For  Case  3,  axial  force  N  is  zero  and  zn  is  constant,  equal  to  0.1981  in. 
**CFS  ■<■  closed-form  solution;  TMS  n,  transfer-matrix  solution. 


Table  6.  Closed-Form  and  Transfer-Matrix  Solutions 
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For  Case  5,  axial  force  N  is  zero  and  zn  is  constant,  equal  to  0.1931 
CFS  "v-  closed-form  solution;  TMS  transfer-matrix  solution. 
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Table  8.  Closed-Form  and  Transfer-Matri x  Solutions 
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For  Case  7,  axial  force  N  is  zero  and  zn  constant,  equal  to  0.1981  in. 
**CFS  ->•  closed-form  solutions;  IMS  ^  transfer-matrix  solutions. 
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closed-form  solutions;  TMS  ^  transfer-matrix  solutions. 
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PART  II 

PREDICTION  OF  BENDING  RUPTURE  STRENGTH 
OF  NON-LINEAR  MATERIALS  WITH  DIFFERENT  BEHAVIOR 
IN  TENSION  AND  COMPRESSION 

Charles  W.  Bert 
The  university  of  Oklahoma 
Norman,  Oklahoma  73019 

Abstract  -  Many  materials  have  quite  different  stress-strain  relations  in 
tension  and  compression.  Examples  include  such  diverse  materials  as  rock, 
cast  iron,  concrete,  tire  cord-rubber,  and  soft  biological  tissues.  It  is 
shown  by  analysis  in  this  paper  that  DBTC  (different  behavior  in  tension  and 
compression)  has  a  profound  effect  on  the  flexural  strength  as  predicted  by 
application  of  fundamental  continuum  mechanics  relations.  The  theory  is 
applied  to  a  non-linear  material  model  which  is  shown  to  be  applicable  to  two 
widely  different  materials:  concrete,  which  has  more  strength  degradation  in 
tension  than  in  compression,  and  steel  cord-rubber,  which  has  strength  which 
is  enhanced  in  tension  by  cord  strengthening  and  degraded  in  -ompression  by 


cord  nicrobuckling 
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1 .  INTRODUCTION 

It  ''as  Ion’  b-.en  known  th.at  certain  materials  behave  significantly 
differeni.lv  in  tension  than  in  compression.  For  example,  this  was  recognized 
by  Saint-Venant  [1],  who  analyzed  the  pure  bending  behavior  of  a  beam  having 
different  non-linear  stress-strain  curves  in  tension  and  compression. 

Specific  experimental  evidence  of  03TC  has  been  reported  for  cast  iron  by 
Gilbert  [21,  cord-rubber  by  Clark  (31,  polymers  by  Zemlyakov  [4],  concrete  by 
Seefried  et  al.  (5],  cortical  bone  by  Simkin  and  Robin  [6],  various  kinds  of 
rock  by  iiaimson  and  Tharp  [7],  and  soft  biological  tissue  by  Pearsall  and 
Roberts  [31.  Of  course,  the  specific  micromechanisns  repons ible  for  D8TC  vary 
from  one  class  of  material  to  another.  Stiff,  brittle  solids  are  weakened  in 
tension  by  microcracks,  while  soft  fiber-reinforced  materials  are  weakened  in 
compression  by  fiber  .nicrobuckling. 

In  this  paper,  a  general  theory  is  presented  for  the  prediction  of  the 
maximum  bending  moment  achievable  in  unidirectional  pure  bending  for  an 
arbitrary  material  with  DBTC .  Subsequently,  the  general  theory  is  applied  to 
a  simple  model  which  is  suitable  for  widely  different  materials:  concrete 
using  a  continuous  damage  mechanics  model  for  the  tensile  behavior,  and 
cord-rubber. 

2.  GENERAL  THEORY 

The  theory  developed  here  is  an  extension  of  the  basic  theory  of 
unidirectional  bending  developed  by  3ach  and  Baumann  (9);  see  Nadai  [10], 
section  22-1.  The  material  is  assorted  to  be  subjected  to  a  unidirectional 
pure  bending  moment,  i.e.,  there  is  no  applied  direct  tension  or  compression. 
The  beam  cross  section  and  loading  are  assumed  to  be  symmetric  about  the  same 
axis,  so  that  bending  takes  place  in  a  single  plane  (no  twisting  or  warping). 
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Poisson-type  lateral  contraction  effects  are  neglected*  and  the  beam  cross 
section  is  assumed  to  be  sufficiently  compact  that  cross-sectional  distortion 
(ovalization  and  warping)  and  buckling  do  not  occur.  Plane  sections  are 
assumed  to  remain  plane  so  that  the  axial  strain  £  is  given 

£  =  <z  ( 1  ) 

where  z  is  the  normal  position  coordinate,  measured  positive  downward  from  the 
neutral-surface  position,  and  <  is  the  bending  curvature. 

The  beam  material  may  have  arbitrary  uniaxial  stress-strain  behavior  in 
the  longitudinal  direction,  i.e., 

<j,.(e)  e  >  0 

o  =  {  (2) 

oc(£)  £  <  0 

where  0  is  the  axial  normal  stress  and  bt(£)  and  oc(e)  are  arbitrary 
continuous  functions. 

In  the  absence  of  direct  axial  loading,  equilibrium  of  forces  in  the 
axial  direction  requires 

0  Ct 

/  ac(<z)b(z)dz  +  /  0t(<z)b(z)dz  =  0  (3) 

-Cc  0 

Here  the  respective  compressive  and  tensile  outer-fiber  distances  are  denoted 
by  Cc  and  Ct  and  the  cross-sectional  width  at  distance  z  is  denoted  by  b(z); 
see  Fig.  1.  Substituting  specific  expressions  for  b(z),  jc(<z),  and  at(<z) 
into  (1),  integrating,  and  using  the  following  geometric  relation  (see  Fig.  1) 
allows  determination  of  Cc  and  Ct: 

Cc  +  Ct  =  H  (4) 

Equilibrium  of  the  internal  moment  due  to  the  stress  distribution  with 
the  externally  applied  bending  moment  M  gives 

‘This  assumption  is  more  reasonable  for  bending  than  for  tension  loading, 
and  more  reasonable  for  brittle  materials,  such  as  concrete  and  rock,  than  for 
ductile  metals  such  as  considered  by  Bert  et  al.  (11). 
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/  zoc  (  <z  )b(  z  )dz  +  /  so  t  ( xz  )  b(  z  )dz  •-  M  (5) 

"Cc  0 

The  ultimate  bending  moment  is  found  by  using  differential  calculus  to 
determine  the  relative  maximum,  i.e. 

M  =  Mmax  when  -lM/dp  -  '  !6) 

where  p  is  any  convenient  parameter  which  can  he  used  as  a  measure  of 
deformation. 

In  principle,  the  general  theory  can  be  applied  to  any  material  {provided 
that  it  has  continuous  but  different*  stress-strain  curves  in  tension  and 
compression)  with  any  singly-symme trie  cross  section.  The  similarity  with  the 
theory  of  ductile  failure  by  plastic  tensile  instability  originated  by 
Considere  [12]  is  readily  .apparent;  c.f.,  Nadai  [10],  section  8-1. 

3.  A  NON-LINEAR  MATERIAL  MODEL 

The  following  mathematical  relationship  is  proposed  to  approximate  the 
stress-strain  behavior  of  a  variety  of  actual  materials: 

a- ( E )  =  E~£  *  ?„e2  ;  e  <  0 

a  =  {  ~  (7) 

at(c)  =  2te  -  F-e2  ;  ►:  >  0 

where  Ec  and  E«.  are  the  Young's  moduli  in  compression  and  tension  and  and 

r*  ^  i  1.3  o  material  constants.  This  model  can  be  considered  to  be  a 

generalization  of  Kragcinovic' s  model  for  plain  concrete  [14].  In  this  case, 


where  is  a  tensi  le-mi crocracx-damage  modulus  (Dr>0)  related  to  .Tanson  and 
Knit's  linear  version  [15]  of  Kachanov's  concept  of  continuum  damage  mechanics 

♦The  case  of  damage  induced  instability  in  beam  bending  for  materials 
having  the  same  continuous  damage  in  tension  and  compression  was  considered 

by  Bos trom  [13]. 
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[16].  For  an  excellent  review  of  CDM,  the  reader  is  referred  to  Hult  [17]. 

The  present  model  is  an  improvement  on  the  Krajcinovic  model  in  that  it 
also  includes  compressive  non-linearity  as  suggested  by  British  Standard  CP1 10 
[18]  : 

Fc  =  -E2/4fcu  (9) 

where  fcu  is  the  material's  ultimate  strength  in  uniaxial  compression  (see 
Appendix  A2). 

To  apply  equation  (7)  to  unidirectional  cord-rubber,  using  experimental 

data  for  cord  rubber  reported  by  Bert  and  Kumar  [19],  one  has  Et  a  Ec,  Fc  =  Ft 

>0. 

If  one  assumes  a  rectangular  beam  of  width  B,  application  of  equation  (7) 
in  equation  (3)  yields: 

J  BK(EtCt2-EcCc2>  +  -j  BK2(FtCt3+FcCc3)  =  0 


or 

3  EcCc2-EtCt2 
K  =  2  FtCtJ  +  FcC~3' 

Integrating  equation  (5),  using  equation  (7),  one  obtains 


(10) 


M  =  J  B<(EtCt3+EcCc3)  +  -j  B<2(FtCt4-FcCc4)  (11) 

Substituting  the  expression  for  <  from  equation  (10)  into  equation  (11) 
and  using  equation  (4)  to  eliminate  Cc,  one  finally  obtains  the  following 
expression  for  M  as  a  function  of  Ct: 


M  =  (3/2) 


[5c(H-Ct)2-EtCt2]  [EtCt3  +  S;;(H-Ct)3] 
FtCt3  + 


[FtCt4  -  Fc  (  H-Ct )  4  ]  [  Ec  (  H-Ct )  2-EtCt2  ]  *" 

+  (  9B/1  6) — — - - - - - - (i2) 

[FtCt3  *  Fc(H-Ct)3]2 


Of  course,  the  ultimate  bending  moment  (Mu)  can  be  found  as  a  function  of 
3 , H, Ec , Et , Fc ,  and  Ft  by  either  direct  substitution  or  use  of  differential 
calculus  (dM/dCt  =  0). 
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4.  APPLICATION  TO  PLAIN  CONCKSTE 

Bert  arui  Xu.na  r  •  201  shoved  that  for  Krajcinovto's  [14]  material  [  Ec  =  Et 
=  »,  Ec  =  0,  and  ?t  given  by  equation  (3)],  equation  (12)  takes  a  form  that 
can  be  expressed  as 

M  =  6  [  1  -  (  1  5/8 )  Kt  *  (3/4).<c2i  1/2)Kt'3]  (13) 

where  ~  o:i/'<OtH"  and  Kt  i  Ct/H.  Sect  and  Kumar  showed  tha t  the  ultimate 
value  of  m,  is  3.407  at  :<c  =  0.5805.  These  values  are  both  slightly 

higher  than  the  values  mu  =  0.354  at  Xt  =  0.5505  obtained  by  Krai ci no vie,  who 
used  the  following  failure  criterion,  which  is  appropriate  only  for  failure  in 
uniform  uniaxial  tension 

da/d£  =0  (14) 

It  is  traditional  to  report  ultimate  strength  data  in  bending  in  terms  of 
bending  rupture  strength  defined  as  the  maximum  tensile  stress  at  ultimate 
bending  moment  calculated  according  to  the  simple,  purely  elastic  formula 

fbu  =  *UC/I  (15> 

where 

I  =  BH3/12  ;  C  =  H/2  (16) 

It  is  *-o  be  emphasised  that  the  value  obtained  by  use  of  equation  (15)  is  not 
a  stress  value;  it  is  just  an  index  of  .no;nent-e  a  crying  ability  expressed  in 
units  of  force  per  unit  area. 

Combining  the  definition  of  mu  (2  5Mu/Dt?W-)  with  equations  (15)  and 
.15),  one  can  eisily  show  that 

^bu  =  ^t^u  (17) 

It  is  shown  in  Appendix  M  that  use  of  th -  criterion  dcr/dc  for  uniform 
axi  il  loading  leads  to  the  following  expression  for  the  ultimate  tensile 
strength: 


(18) 
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Thus, 

^bu/“  tu  =  4t[tu  (19) 

Fig.  2  shows  a  comparison  among  the  experimental  data  of  Gonnerman  and 
Shuman  [21]  and  of  Brooks  and  Neville  [22]  (both  wet  and  dry  stored)  and  the 
theoretical  prediction  of  Krajcinovic  [14]  and  of  Bert  and  Kumar  [201. 

It  is  noted  that  the  Krajcinovic  and  Bert  and  Kumar  predictions  both  are 
fbU  vs.  ftu  relations  that  are  straight  lines  starting  at  the  origin  (with 
slopes  of  1.42  and  1.63,  respectively).  However,  the  experimental  curves  are 
straight  lines  shifted  upward  so  that  they  have  f^u  intercepts.  Apparently 
this  upward  shift  is  due  to  statistical  considerations,  such  as  considered  by 
Weil  and  Daniel  [23],  for  instance,  using  the  Weibull  distribution  [24]  and  by 
Jayatilaka  [25],  using  a  different  approach. 

It  is  well  known  that  concrete  has  a  stress-strain  non-linearity  in 
compression,  although  less  pronounced  than  the  one  in  tension.  This  was 
recognized  in  British  Standard  CP110  [18];  see  also  Refs.  [26-30].  This 
non-linearity  is  represented  by  the  term  in  equation  (7)  with  a  coefficient 
Fc.  To  obtain  a  quantitative  value  for  Fc,  one  must  know  either 

(1)  the  mean  compressive  stress-strain  curve  to  failure 

(2)  the  mean  compressive  strength  (fcu),  for  it  is  shown  in 
Appendix  A2  that  Fc  =  E^/Dc  where  Dc  =  4fcu. 

(3)  the  rati;-  of  mean  compressive  strength  to  mean  tensile  strength 
for  Dc/Dc  =  fcu/ftu* 

Since  the  author  has  not  yet  been  able  to  obtain  a  copy  of  Ref.  [21],  it 
is  necessary  to  consider  typical  values  of  fcu/ftu*  (Sef.  [22]  had  fcu/ftu  “ 
15.)  From  statistical  considerations,  for  various  values  of  the  Weibull 
parameter  and  sample  size,  Jayatilaka  [25],  Table  5.13,  gave  values  of  fcu/ftu 
ranging  from  3  to  20.  However,  the  ratio  fbu/^tu  predicted  by  the  present 
theory  is  not  strongly  affected  by  fCu/*tu»  as  can  be  seen  in  Table  1. 
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Piechnik  and  Pachla  [31]  proposed  a  power-law  damage  function  to  describe 
the  tensile  behavior  of  plain  concrete.  An  analysis  using  this  damage 
function  to  predict  failure  in  bending  an!  in  tension  is  presented  in  Appendix 
B.  For  n  -  4,  the  value  recommended  by  Piechnik  and  Pachla,  Uv^/D  =  mu  = 

0.715  at  a  value  of  0.529  for  Kt  foe  the  bending  case.  For  tension,  a  value 
of  0.669  is  attained  for  otu/ D.  Thus, 

Jbu/°tu  =  5/0.669  -  1.07 

This  value  of  1.07  is  considerably  lower  than  is  consistent  with  the 
experimental  results  of  Refs.  [21]  and  [22].  The  explanation  for  this 
discrepancy  may  be  any  one  of  the  following  factors  or  a  combination  of  them: 

(1)  Difference  in  details  of  cement,  composition  and  fineness, 
aggregraie  sice,  material  quality,  and  curing  [32]  between  the 
1925  American  concrete  [21]  and  the  1977  British  concrete  [22] 
on  one  ha  no  and  the  1979  Polish  concrete  [31]  on  the  other. 

(2)  Inert  se  in  measured  ^bu/^tu  ^ue  to  a  brittle  failure 
mechanism  and  statistical  sice  effects  as  discussed 

r_  r  •• ,  icusiy . 

(3)  Possibility  that  the  material  damage  relation  for  concrete 
has  a  threshold  effect,  as  suggested  by  Young  [33]  and  Bert 
and  Kumar  [20],  rather  than  the  gh  power-law  damage  function 
suggested  in  [31], 

5.  APPLICATION  TO  CORD- REINFORCED  RUBBER 

It  was  shown  m  the  recent  experiments  of  3ert  and  Kumar  [19]  that  the 
stress-strain  relations  of  this  class  of  materials  is  strongly  dependent  upon 
the  type  of  cord.  For  example,  arami 1-rubber  loaded  in  the  cord  direction 
(9°)  has  i  irastic  discontinuity  in  slope  at  the  origin  ( E~  >  Ec ) ,  yet  the 
same  material  loaded  at  JO*  to  the  cord  direction  is  very  nearly  linearly 
through  ( Ec  =  Et,  Fc  =  ?►.  =  9).  As  an  example  of  the  application  of  the 


i 
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generalized  model  presented  in  Section  3,  steel  cord-rubber  specially  made  to 
have  highly  waved  cords  is  considered.  Then  Ec  “  Et  =  E  and  Fc  =  Ft  =  F  >  0. 
Then  equation  (12)  simplifies  to 

2 

M/BH2  =  (  E2/F)  (  1  -2Kt )  [  3  -  |^(  1  -2Kt+2Kt2  )  ( r— ^ )  ]  (20) 

3K^-3Kt+1 

It  can  be  shown  that  the  maximum  vlaue  of  dimensionless  bending  moment  is 
6MUF/BE2H2  =  0.433  at  a  value  of  0.380  for  Kt.  Now  due  to  the  monotonically 
increasing  slope  of  the  tension  portion  of  the  stress-strain  curve,  it  does 
not  exhibit  a  maximun-load  phenomenon.  (Tensile  failure  is  obviously  due  to 
another  mechanism,  probably  cord  failure.)  Thus,  in  this  case,  it  is 
necessary  to  ratio  the  bending  strength  to  the  compressive  strength  rather 


than  the  tensile  strength.  From  Appendix  A2,  one  has 

acu  =  -E2/4F 

Thus, 


(21  ) 


6MU/BH2 

latu/°cu!  ~  2/1-  =  1,73 

EV4? 

This  indicates  that  abu  is  1.73  times  the  compressive  strength.  In  reality, 
the  tension  portion  of  the  steel-cord  rubber  is  even  steeper  than  indicated  by 
a  power  of  two.  Thus,  abu  would  be  expected  to  be  even  greater.  This  could 
not  be  verified  experimentally  due  to  the  extremely  flexible  nature  of  the 
matieral,  i.e.  its  flexural  stiffness  is  too  low  to  measure  in  practice. 


6.  CONCLUSIONS 

A  general  theory  was  developed  for  predicting  the  maximum  bending  moment 
that  can  be  carried  by  beams  having  arbitrary  compact  cross  sections  and 
constructed  of  material  havinq  an  arbitrary  non-linear  stress-strain  curve 
that  is  different  in  tension  and  in  compression.  As  an  application  of  the 
theory,  it  was  applied  to  rectangular-cross-section  beams  of  plain  concrete 


and  of  cord- rubber 
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A? PEND XX  A:  RELATIONS  BETWEEN  TENSILE  AMD  COMPRESSIVE 
STRENGTHS  AND  MATERIAL  DAMAGE 

1 .  Tension  Loading 
Using 

0  =  E£  +  Ft£2  ;  £  0  0 

it  is  obvious  that  material  instability  cannot  occur  unless  ?t  <  0. 

v 

Considere  criterion  [12]  is 

do 


de 


=  0 


(A!  ) 

Then  the 

( A2 ) 


or 


E  r  2FtC  ~  0 

rhus,  the  strain  at  which  tensile  instability  occurs  is 


-ti 


-E/2F 


and  the  tensile  strength  is 


(A3) 


'tu 


2Ft  4Ft 


4F 


=  Dt/4 


( A4 ) 


2.  Compression  Loading 
Using 

a  =  Se  +  Fc£2  ;  £  <  0  (A5) 

it  is  obvious  that  material  instability  cannot  occur  in  this  situation  unless 
Fc  >  0.  Again  using  the  Considere  criterion,  equation  (A2),  or.e  obtains  the 
following  expression  for  the  i.nstabilty  strain: 

sci  =  -  E/2F,  ( A6 ) 


and  the  predicted  compressive  strength  is 


64 


o 


cu 


S2/4Fc  =  Dc/4 


(A7) 


APPENDIX  B:  EFFECT  OF  POWER-LAW  TENSILE- DAMAGE 
FUNCTION  ON  BENDING  STRENGTH  OF  CONCRETE 


1 .  Bending  Loading 

In  the  continuum  damage  theory  [17],  the  nominal  stress  a  is  related  to 
the  net  or  true  stress  3  by  the  following  relation 

0  =  (1-u>)S  (31  ) 

Here  u  is  called  the  damage. 

Piechnik  and  Pachla  [30]  proposed  the  following  power-law  damage  function 
for  concrete  in  tension: 

ui  =  (S/Dt)n  (B2) 

Here  n  and  Dt  are  material  constants.  Based  on  limited  test  data,  Piechnik 
and  Pachla  recommended  a  value  of  4  for  n.  Thus,  the  following  expressions 
for  to  are  used  here: 


(S/Dt)n  for  S  >  0 

W  =  {  (B3) 

0  for  3  <  0 

Since 


3  =  EE  ( B4 ) 

we  obtain  the  following  expression  for  the  stress  distribution 

<JC(E)  =  EE  E  <  0 

a  =  {  (35) 

at(E>  =  EE  -  E(E/Dt)nen+1  £  >  0 

It  is  readily  apparent  that  this  model  is  a  generalization  of  the 
linear-damage  model  hypothesized  by  Krajcinovic  [14]. 

Use  of  equation  (B5)  in  equation  (3)  leads  to  the  following  relationship 
between  the  curvature  <  and  the  dimensionless  tensile-outer-fiber  distance  (Kt 
=  Ct/H): 


<  = 


(Dt/EHKt)  [1+(n/2)  ]  1/r'(2Kc-1  -  Kt-2)1/" 


( B6 ) 
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Similarly,  substitution  of  equation  (35)  into  equation  (5)  and 
integration  gives 

3M/3SKH3  =  Kt3+(1-Kc)3  -  !3(2+n)/2(3+n) ] (2KZ2-Kt)  (37) 

Finally,  substitution  of  equation  (36)  into  equation  (37)  yields 
m  =  6M/BOcH2  -  2  [  1  +( n/2  )  ]  1  -''n  (  2Kt"  !  -<t-2  )  1  -/n  [  Kt~  '  -  3+N  +  (3-2N)K.l  (36) 

Obviously,  for  the  case  of  linear  damage  (n  =  1),  equation  (33)  reduces  to 
equati on  (13). 

2.  Tension  Loading 

Applying  the  tensile  instability  criterion,  equation  (A2)  to  ot(S)  as 
given  by  equation  (85),  one  obtains  the  following  expressions  for  the 
instability  strain  and  ultimate  tensile  strength: 

eti  =  (Dt/2)(n+l)_,/n  (39) 

0tu  =  D(n+n-1/n[i  _(n  +  i  ) -ntl  ]  (310) 

Table  1.  Effect  of  ratio  of  compressive  to  tensile  strength  on 
ratio  of  bending  to  tensile  strength 


f"  /  f 

CCU’  -  tu 


1  0 
3 


fbu/ftu 


0 

1  .63 

1/10 

1 .61 

1/3 

1.57 
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